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ABSTRACT 

A large population of Intermediate Mass Black Holes (IMBHs) might be produced 
at early cosmic times as a left over of the evolution of the very massive first stars. 
Accretion onto IMBHs provides a source of (re)ionizing radiation. We show that the 
baryon mass fraction locked into IMBHs and their growth is strongly constrained by 
the observed residual Soft X-ray Background (SXRB) intensity. Thus, unless they are 
extremely X-ray quiet, miniquasars must be quite rare and/or have a short shining 
phase. As a byproduct, we find that miniquasars can not be the only source of reion- 
ization and that their alleged contribution to the near infrared bands is completely 
negligible. 

Key words: black hole physics - galaxies: formation - intcrgalactic medium - cos- 
mology: theory - diffuse radiation 



1 INTRODUCTION 

Recent numerical (Bromm, Coppi & Larson 1999, 2002; 
Abel, Bryan & Norman 2000, 2002) and semi-analytical 
(Omukai & Nishi 1998; Omukai 2002; Nakamura & 
Umemura 2001, 2002; Schneider et al. 2002, 2003; Omukai 
& Palla 2003) studies consistently predict that the first, 
so-called Population III (PopIII), stars had characteristic 
masses of 100-600 Mq, i.e. roughly ten times more mas- 
sive than those observed today. Since mass loss at zero 
metallicity can be neglected (Kudritzki 2002), the final fate 
of PopIII stars is essentially set by their initial mass. In 
the narrow range of progenitor masses, 140 Mq < Mi, < 
260 Mq , PopIII stars explode as pair instability supernovae 
(SN-y^y) leaving no remnant and polluting the universe with 
metals (Heger & Woosley 2002). Non-rotating stars with 
masses < 140 Mq and > 260 Mq will instead collapse di- 
rectly into black holes (BH), without metal ejection. These 
intermediate-mass BHs (IMBHs) may provide the seed for 
the supermassive BHs observed in the center of luminous 
galaxies (e.g., Volonteri, Haardt & Madau 2003). If IMBHs 
accrete material from the surrounding medium, they will 
shine as "miniquasars" and may contribute to the reioniza- 
tion of the Universe at high redshift (Haiman & Loeb 1998; 
Venkatesan, Giroux & ShuU 2001; Oh 2001; Ricotti & Os- 
triker 2004a,b; Madau et al. 2004). 

In this paper, we discuss the contribution of this popula- 
tion of early miniquasars to various cosmic backgrounds. We 
show that available observations can be used to set strong 
constraints on the radiative proprieties of these sources. In 
particular the Soft X-ray Background (SXRB) is easily over- 
produced unless miniquasars are very rare, or, in case more 



common, if their shining phase is very short. As a byprod- 
uct, we show that, in both cases, their expected contribution 
to the Near Infrared Background is negligible. 



2 COSMIC BACKGROUND RADIATION 

The mean specific background intensity J{uo, zo) at redshift 
zo observed at frequency fo, produced by a population of 
sources characterized by a comoving emissivity jiy{z), can 
be written as 



(l + ^o)' 
4-K 



(1) 



where u — i'o{i + z) / (1 + zq) , and dl/dz is the proper line ele- 
ment. In the above expression we have neglected, for reasons 
that will become clear later on, any absorption of photons as 
they propagate in the expanding Universe. The source term 
ji, can be written, in general, at any given cosmic time t, as 
a convolution of the light curve of the sources lv{t) with the 
source formation rate (per unit comoving volume) dN/dt: 



dN 



Ut^t')^dt' ^Irsie^'f/^^ 



,dN , ^ 

D^W. (2) 



where rs = eMqc / JeLe ^ 4.4 x 10 e/fs yr is the 
Salpeter time, e is the accretion radiative efficiency, and 
Le 2i 1.3 X 10^**(A/i3H/MQ) is the Eddington luminosity. 
The second approximated equality holds averaging the light 
curve over a typical source lifetime r;/, leading to a spe- 
cific luminosity Z^. and assuming a constant source forma- 
tion rate over such timescale. The emissivity is then pro- 
portional to the accreted mass AMacc = (e^'^''^^ — 1)Mbh, 
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where Mbh is the initial BH mass. Although the frequency 
integrated background depends only on the amount of ac- 
creted mass, and not on the details of the accretion process, 
the spectral energy distribution (SED) in a given narrow 
band may depend on the properties of the accreting BHs. 
For example, for fixed AMacc, a shorter t;/ would imply a 
smaller final BH mass (and a larger number of those BHs) , 
with a resulting hotter disk component (see next Section 3) . 

Note that, since can be written as a fraction /b of 
the Eddington luminosity (see next Section), ./(z/o, zo) is in- 
dependent of this parameter. We compute the source for- 
mation rate using the Press-Schechter formalism (Press & 
Schechtcr 1974) adopting the minimum dark matter halo 
mass, Mmin[z), computed by Fuller & Couchman (2000). 

Throughout the paper, we adopt the 'concordance' 
model values for the cosmological parameters: h = 0.7, 
Qm = 0.3, Qa = 0.7, ^6 = 0.044, trg = 0.9, and F = 0.21, 
where h is the dimensionless Hubble constant. Ho = lOO/i 
km s"^ Mpc~^ 



3 MINIQUASAR SPECTRUM 

The physical characterization of the sources is encoded in 
the SED Iv, and in the typical lifetime r. The UV/X-ray 
SED of an accreting BH, is, observationally, approxima- 
tively described in terms of two main continuum compo- 
nents (see, e.g., Tanaka & Levin 1995). The low energy com- 
ponent is thought to originate from the putative accretion 
disk, and, at least in stellar-sized BHs (the so-called galac- 
tic black hole candidates GBHCs), it is spectroscopically 
well described by a "multicolor disk black body" (Mitsuda 
et al. 1984). The accreting gas is assumed to be optically 
thick, and the gravitational power locally released as black 
body radiation. Different radii radiate at different tempera- 
tures, with the hottest Planckian emitted at ~ 5i?s, where 
Rs = 2GM/c? . Assuming Eddington limited accretion, this 
model yields to kT^ax « 1 hiV{MBH /Mq)'^'^ (Shakura & 
Sunyaev 1973). The characteristic multicolor disk spectrum 
is broadly peaked at Speak ^ ^kTmax, follows a power law 
with lv,MCD oc v^f'^ for energies hv < Speak, and exponen- 
tially rolls off for E > Speak- 

The second observed component is a simple power-law 
lv,PL oc with a ~ 1. The precise origin of this power- 
law emission is controversial (several difi^erent models exist 
to date), and is likely to originate from synchrotron and/or 
inverse Compton (IC) emission by a mixture of thermal and 
nonthermal electrons, located in an active corona and/or in a 
relativistic jet (for e recent review see, e.g., Zdziarski & Gier- 
linski 2004). Both the low and the high energy ends of the 
power law are not firmly known. Hard-states of GBHCs and 
Seyfert I galaxies show exponential cut-offs at S ~ 100 — 500 
keV (Matt 2001), while in the steeper soft state of GBHCs, 
and in more distant objects such as QSOs, observations fail 
due to the low S/N. The low energy end of the power law, 
as well, could start at E ~ Speak, if the soft photon in- 
put for the IC cooling of the relativistic electrons is due to 
disk radiation, or, instead, could extend into the IR if lo- 
cal synchrotron radiation is, in the corona/jet, energetically 
important (SSC models). 

The two component spectrum, though historically mo- 
tivated by studies of GBHCs and Seyfert galaxies, approx- 



imates quite well the high energy emission observed in ul- 
traluminous X-ray sources, a recently discovered population 
of X-ray sources possibly associated to accreting IMBHs, 
Mbh ^ 500 - 10^ M© (Miller et al. 2003). Observation- 
ally, an IMBH seems to have comparable MCD and PL lu- 
minosities, although the ratio between the two components 
might vary among sources, and the extension of the PL is 
not known. In order to model in a simple, but physically 
motivated, way the SED of an accreting IMBH, we assume 
that the PL component is due to IC scattering of thermal 
disk photons, and hence we consider energies S > Speak- 
This is a clear difference with respect to the SED used by 
Dijkstra ct al. (2004), and by Madau ct al. (2004), who as- 
sumed the PL component to be present for S > 13.6 eV, 
so that it is dominant with respect to the MCD in terms of 
emitted ionizing photons. 

We account for the population spread of the PL/MCD 
flux ratio by introducing an empirical parameter: 

LmCD J li,,MCDdl' 

It may worth noticing that, for a fixed BH mass, the value of 
$ determines the luminosity of the X-ray part of the emitted 
spectrum. We also parameterize the bolometric luminosity 
(i.e., MCD-I-PL) of an IMBH of mass Mbh as a fraction fs 
of the Eddington luminosity, Le ^ 1-3 x W^^{Mbh/Mq) 
erg s~^- 



4 THE SOFT X-RAY BACKGROUND 

Moretti et al. (2003) determined the intensity of the total 
SXRB in the energy range 0.5-2 keV to be (7.53 ± 0.35) x 
10~^^ erg cm~^ deg~^ when combining 10 different mear 
surements reported in the literature. Including further deep 
pencil beam surveys together with wide field shallow sur- 
veys, they find that 94l^% of the SXRB is made up of dis- 
crete X-ray sources (the majority being point sources) at low 
redshift, z < A (Barger et al. 2002, 2003). By reanalyzing the 
Moretti et al. uncertainty budget, Dijkstra, Haiman & Loeb 
(2004) provided a mean and a maximum intensities of the 
unaccounted SXRB flux, 0.35 x 10"^^ and 1.23 x 10"^^ erg 
cm~^ s~^ deg~^, respectively. 

A population of IMBHs forming at high redshift can 
contribute to the SXRB. We compute here the expected 
background intensity in the energy band 0.5-2 keV accord- 
ing to Eqs. (l)-(2) for a population of miniquasaxs whose 
formation continues down to a given redshift, Zend, as a fi- 
nal evolutionary product of massive PopIII stars. Given the 
hard energy band and low IGM metallicities we are con- 
cerned with, we neglect any absorption term in the radiative 
transfer equation (eq. 1). 

In the 0.5-2 keV energy range, the background intensity 
is dominated by the power-law component of the miniquasar 
spectrum, unless "1> <C 1. We consider here the two extreme 
values of $ = 1,10~^. In order to set upper limits on the 
propriety of this first population of accreting IMBHs, we 
give the result for the maximum unaccounted SXRB flux 
derived by Dijkstra et al. (2004) . 

Unless miniquasars are extremely X-ray quiet (and 
therefore $ <C 1), we find that the SXRB sets strong con- 
straints on the density of miniquasars, which are summar- 
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Figure 1. Constraints on the IMBH lifetime, th, and fliMBH /^b 
for redshift ^e^d = 6 {solid line), z^^^ = 12 (dotted), z^^d = 18 
(dashed), and z^^^ = 24 (long- dashed). The top (bottom) panel 
refers to the case <I> = Lpi^/ Lmcd = !('!• = 10~^). We adopt 
here e = 0.1. 
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Figure 2. Maximal growth of IMBHs as function of the initial 
mass density of IMBHs forming at 2 = 24. Different lines refer 
to different value of e. Top panel: final density in BHs. Bottom 
panel: growth factor. 



rizod in the top panel Fig. 1 as a function of their life- 
time. The curves refer to different turn-off redshifts z^nd = 
6, 12, 18, 24, respectively. Apart from the differences in the 
value of Zend, which introduces an uncertainty of a factor w 2 
on the estimates, we see that for $ = 1, the mass fraction of 
IMBH cannot exceed lO^^f^b (that is of the same order of 
the density of the SMBH today, Merritt & Fcrrarese 2001), 
even for an extremely short lifetime tu — 10^ yr. This con- 
straint increases to 0.imbh < 0.1 ilf, in the case in which 
IMBHs are extremely ineflBcient X-ray emitters ($ — 10~^), 
i.e. a large fraction of the baryon density might be locked 
into BHs without exceeding the SXRB constraint (bottom 
panel of Fig. 1). Since we have assumed the maximal SXRB 
residual intensity, the above values must be understood as 
strong upper limits. In addition, other sources, such as high 
redshift quasars (Dijkstra ct al. 2004), may contribute to the 
unresolved SXRB, leading to even more stringent limits . We 
conclude that early miniquasars were quite rare and/or their 
shining phase lasted only for a very short period of time. 

Wc can rewrite those limits in terms of the maximal 
mass growth of BHs allowed by the unaccounted SXRB. 
In Fig. 2 we show the final IMBH density and the growth 
factor, defined as the ratio between final and initial IMBH 
density, for sources forming at 2; = 24 (different lines refer to 
different values of accretion radiative efficiency e) as function 
of the initial mass density of IMBH, QiMBH.in- For sources 
forming down to Zend = 12 {zend = 6) the limits are tighter 
of a factor 1.8 (2.6). 

We find that for e = 0.1 a strong upper limit to the 
final mass density 7 x 10^^ for a wide range of initial 
densities, i.e. for ^Iimbh /^b ^ 10~® Considering that 
low redshift {z < 6) accreting BHs are taken into account in 
the resolved fraction of the SXRB, we can derive an upper 



limit on the density of BHs (active and not) at 2; > 6 of 
piMBH < 3.8 X 10"' Mq Mpc"^, which is ~ 10% of present 
day SMBH mass density (Yu & Tremaine 2002). This value 
is not at odd with current models of SMBH assembly in 
a hierarchical structure formation scenario (Volonteri et al. 
2003; Madau ct al. 2004). As example, Volonteri et al. (2003) 
find that the mass density locked into BHs is of the order 
of 10'' Mq Mpc~^ aX z ^ 6. This mass growth corresponds 
to a SXRB contribution of 0.29 x 10"^^ erg cm"^ s"^ 
dog^^, that is a not negligible fraction of the maximal unre- 
solved intensity, and it is comparable to the mean value of 
this quantity. Moreover, we find also a contribution in the 
hard-X band (2-10 keV, so-called HXRB) of the order of 
0.33 X 10~^^ erg cm~^ s"'^ deg~^, corresponding to ~ 1.7% 
of the observed HXRB. Thus, within hierarchical clustering 
models, a significant fraction of the unaccounted SXRB (and 
HXRB) should come from the growth of IMBHs in the early 
Universe. For example, using the Madau et al. (2004) model 
we find that the SXRB (HXRB) is ~ 0.20 (0.23) x 10"^^ erg 
cm~^ s~^ deg~^ already at 2; = 14 for 3.5a peak seeds, as- 
suming an accreted mass corresponding to 10~^ of the halo 
mass. 

Stronger constraints on the maximal growth of IMBHs 
in the early Universe can be set by the next generation of X- 
ray satellites (e.g. Constellation-X, XEUS) that will be able 
to resolve sources 10 times fainter than the present facilities. 
Extrapolating the log N/ log S to this flux limit will allow to 
resolve the SXRB entirely (Moretti et al. 2003). If this were 
the case, strong limits on the grow history of BHs in the 
early Universe can be derived, or conversely, on the spectral 
energy distribution at high energies of these sources. 
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5 LIMITS ON REIONIZATION 

The limits on the radiative proprieties of miniquasars de- 
rived from the SXRB constraint allow us to give a simple 
estimate of the role of these sources in the reionization of 
the Universe. 

The number of ionizing photon per hydrogen atom can 

be written as 

nion _ ^TUhC'^ fuV ( ^IMBh\ /.n 
■^~^^(MV 0.b J fir, ^' 

where {^imbh /^b) fin = {e'^'f^'^^ - 1){CIimbh /i^b)i'n, fuv 
is the fraction of the bolometric power emitted as hydrogen- 
ionizing photons with mean energy (hi/) , mn is the hydrogen 
mass, and X ~ 0.76 is the mass fraction in hydrogen. In Fig. 
3 arc shown the limits on this quantity derived from the 
maximum unaccounted SXRB assuming e = 0,1. Different 
lines refer to different lifetime of the miniquasar phase, t;/. 
The labels report the number of e- folding times, r;//rs. 

Sources forming at very high rodshift and accreting for 
10 e-folding time will be able to produce just 3 photons per 
hydrogen atoms, but shorter lifetimes give lower nion/nn, 
indicating that miniquasars can not easly reionize the Uni- 
verse if recombination is taken into account. At lower red- 
shift the situation is even worst. Miniquasars forming down 
to redshift 9 can not produce more than one ionizing photon 
per hydrogen atom without saturating the SXRB. Moreover, 
these values are to be taken as strong upper limits, since we 
considered the maximum residual SXRB intensity. Using the 
mean unaccounted SXRB intensity will lead to a reduction 
of a factor ~ 1/3 of rnon/nH- 

In conclusion, miniquasars unlikely account for the 
reionization of the Universe even at high rodshift, unless they 
are extremely X-ray quiet. In order to have nion/TiH ~ 10 
without exceeding the maximum unresolved SXRB we must 
required $ < 0.15 (0.07), for sources forming down to red- 
shift z = 24 {z = 9) and living t;/ ~ 4ts . 



6 THE NEAR INFRARED BACKGROUND 

Recent measures of the Near Infrared Background (NIRB; 
sec Hauscr & Dwck 2001 for a review) show an excess 
with respect of the observed light from galaxies in deep 
field surveys (Madau & Pozzetti 2000; Totani et al. 2001). 
The discrepancy is maximal in the J band corresponding to 
1.7 — 4.8 X 10~* erg s"'^ cm~^ sr~^. The large uncertainty 
on this value is given by the different adopted subtraction of 
zodiacal light (i.e. sunlight scattered by the interplanetary 
dust) contribution*. 

Estimates based on theoretical models suggest that this 
excess can be well produced by redshifted light from the first 
very massive (M > 100 M©) PopIII stars (Santos, Bromm 
& Kamionkowski 2002; Salvaterra & Ferrara 2003) if these 
form efficiently down to Zend = 9. The same stars can also 
account for the observed small scale angular fluctuations 
detected in the same bands (Magliocchetti et al. 2003). In 
order to avoid over-enrichment of the IGM with metals at 

* The lower limit is obtained for the zodiacal light model of 
Wright (1998), whereas the upper limit is for the Kelsall et al. 
(1998) one. 
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Figure 3. Limits on the number of ionizing photons per hydro- 
gen atom as function of redshift given by the unresolved SXRB. 
Different lines refer to different lifetime, Tif. The labels report the 
number of e-folding Tif/rs- We adopt here e = 0.1. 

high redshift, most of these massive stars must end up into 
IMBHs. Cooray & Yoshida (2004) have speculated that if 
these IMBHs accrete matter and shine as miniquasars, they 
might give an important contribution to the NIRB. 

We have revisited this conclusion in the light of the 
results of the Section 4. Using the limits implied by the 
SXRB, we find that the contribution to the NIRB of these 
sources is completely negligible. In fact, for Zend = 9, the 
NIRB contribution from miniquasars in the ,J band allowed 
by the SXRB constrain is ^ 10~® erg cm~^ sr~^, hence 
well below the observed value. 

An appreciable contribution in the NIR bands is pos- 
sible only in the unlikely case in which miniquasars are ex- 
tremely X-ray quiet. In this case, $ <g 1 and the limits set 
by the SXRB are very weak. For sources forming down to 
Zend = 9 the NIRB excess data can be fitted without satu- 
rating the unaccounted SXRB. By redshift 9 and considering 
e = 0.15, BHs can increase their mass density for almost 6 
e-folding times, or about a factor 400. In this case, an initial 
BH mass density of 10~* Q,b is sufficient to reproduce the 
NIRB data. On the other hand these BH have to accrete all 
the time down to 2 = 9 resulting in a final mass density of 
niuBH ~ 0.07 rib. 



7 DISCUSSION 

We have studied the contribution of the first generation of 
miniquasar to cosmic backgrounds. In particular, we have 
shown that the observed residual SXRB intensity (Moretti et 
al. 2003, Dijkstra et al. 2004) can be used to set strong con- 
straints on the abundance and radiative efficiency of these 
sources. Unless these objects are extremely X-ray quiet, the 
SXRB is easily overproduced, requiring miniquasars to be 
quite rare and/or have a short shining phase. Should accret- 
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ing IMBHs saturate the SXRB, they would contribute also 
6-7% of the HXRB. 

As a consequence of our analysis, it is unlikely that 
miniquasar can reionize the Universe, since their are limited 
to produce < 3 photons per hydrogen atom, oven at high 
redshift. This conclusion is similar to Dijkstra et al. (2004), 
though our limits are tighter owing to a more physically 
motivated description of the miniquasar spectrum. More- 
over our approach allow us to follow the evolution of the 
mass density of accreting IMBHs with time, so that we can 
derive important constraints on the mass growth of these 
objects in the early Universe. We derived a strong upper 
limit to (active and not) IMBHs mass density at z > 6, be- 
ing pBH < 3.8 X 10"' Mq Mpc"'\ Although this value is not 
at odd with current model of SMBH assembly in the hierar- 
chical scenario of structure formation, stronger constraints 
on the SXRB unaccounted fraction by future X-ray facili- 
ties (i.e. ConstoUation-X, XEUS) could question our ideas 
of the formation of quasars. In fact, given the prediction of 
these models we expect that a not negligible fraction of the 
SXRB will not be resolved, being the signature of the growth 
of IMBHs in the early Universe. 

As a furtlu-r byproduct, wo have shown that their contri- 
bution in the near infrared bands is completely negligible. 
In the proposed models of the NIRB (Santos et al. 2002; Sal- 
vaterra & Ferrara 2003, Maghocchetti et al. 2003) the NIRB 
excess is due to redshifted light of PopIII stars with masses 
larger than 100 Mq. In order to avoid over-enrichment of 
the IGM at high redshift, most of these stars must end up 
in IMBHs, locking « 10% of the baryons into compact ob- 
jects already a,t z — 9. Though not excluded by any of the 
current experiments (including gravitational lensing data, 
Wambsganns 2002), this requirement is somewhat extreme, 
as pointed out by Madau & Silk (2005). On the other hand, 
we have shown that miniquasars powered by accretion onto 
IMBH cannot contribute appreciably to the NIRB, as they 
easily exceed the SXRB constraints. In fact, the IMBHs 
left over by the first stars must be characterized by a very 
short shining phase (< 10^/^^ yr, assuming $ « 1) and/or 
very low accretion efficiency in order not to overproduce the 
SXRB. As a consequence, IMBHs carmot grow appreciably 
in mass. Only in the unrealistic case $ <C 1, the contribution 
to the unaccounted NIRB from miniquasars might dominate 
that of the progenitors. In this case, we found that ~ 7% of 
the baryons must be locked into IMBHs at 2 = 9. 
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